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Abstract 

The notions of spectral stability and the spectrum for the Vlasov- 
Q\ Poisson system linearized about homogeneous equilibria, fo(v), are 

reviewed. Structural stability is reviewed and applied to perturbations 
of the linearized Vlasov operator through perturbations of /o- We 
prove that for each fo there is an arbitrarily small 5/q in W ' (K) such 
O that /o + 5fo is unstable. When fo is perturbed by an area preserving 

rearrangement, fo will always be stable if the continuous spectrum 
T"! is only of positive signature, where the signature of the continuous 

, ^ spectrum is defined as in [H [2] . If there is a signature change, then 

^ there is a rearrangement of fo that is unstable and arbitrarily close to 

fo with /q in W 1 ' . This result is analogous to Krein's theorem for the 
continuous spectrum. We prove that if a discrete mode embedded in 
the continuous spectrum is surrounded by the opposite signature there 
is an infinitesimal perturbation in C n norm that makes fo unstable. 
If fo is stable we prove that the signature of every discrete mode is 
the opposite of the continuum surrounding it. 
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1 Introduction 



The perturbation of point spectra for classical vibration and quantum me- 
chanical problems has a long history [31 Hj. The more difficult problem of 
assessing the structural stability of the continuous spectrum in e.g. scatter- 
ing problems has also been widely investigated E]. Because general linear 
Hamiltonian systems are not governed by Hermitian or symmetric operators, 
the spectrum need not be stable and a transition to instability is possible. 
For finite degree-of-freedom Hamiltonian systems, the situation is described 
by Krein's theorem El E], which states that a necessary condition for a 
bifurcation to instability under perturbation is to have a collision between 
eigenvalues of opposite signature. The purpose of the present paper is to 
investigate Krein-like phenomena in Hamiltonian systems with a continuous 
spectrum. Of interest are systems that describe continuous media which are 
Hamiltonian in terms of noncanonical Poisson brackets (see e.g. [TUI ITTj ). 

Our study differs from that of [T2], which considered canonical Hamil- 
tonian systems with continuous spectra in a Hilbert space where the time 
evolution operator is self-adjoint. The effects of relatively compact pertur- 
bations on such a system were studied and it was proved that the existence 
of a negative energy mode in the continuous spectrum caused the system to 
be structurally unstable. It was also proved that such systems are otherwise 
structurally stable. Also, our study differs from analyses of fluid theories 
concerning point spectra [131 E] an d point and continuous spectra [T^j, the 
latter using hyperfunction theory. 

A representative example of the kind of Hamiltonian system of interest is 
the Vlasov- Poisson equation [IB], which when linearized about stable homo- 
geneous equilibria gives rise to a linear Hamiltonian system with pure con- 
tinuous spectra that can be brought into action-angle form [H [TTJ, HH1 E] ■ A 
definition of signature was given in these works for the continuous spectrum. 
In the present paper we concentrate on the Vlasov-Poisson equation, but 
the same structure is possessed by Euler's equation for the two-dimensional 
fluid, where signature for shear flow continuous spectra was defined [191 EQ] 
and, indeed, a large class of systems [21J. Thus, modulo technicalities, the 
behavior treated here is expected to cover a large class of systems. 

In Sec. [2] we review on a formal level the noncanonical Hamiltonian struc- 
ture for a class of systems that includes the Vlasov-Poisson equation as a 
special case. Linearization about equilibria is described, the concept of dy- 
namical accessibility, and the linear Hamiltonian operator T, the main sub- 
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ject of the remainder of the paper, are defined. In the remainder of the paper 
we sketch proofs, in varying levels of detail, pertaining to properties of this 
linear operator for various equilibria. In Sec. [3] we describe spectral stability 
in general terms and analyze the spectrum of T for the Vlasov case. The 
existence of a continuous component to the spectrum is demonstrated and 
Penrose plots are used to describe the point component. In Sec. [4] we describe 
structural stability and, in particular, consider the structural stability of T 
under perturbation of the equilibrium state. We show that any equilibrium is 
unstable under perturbation of an arbitrarily small function in W 1,1 . In Sec. [5] 
we introduce the Krein-Moser theorem and restrict to dynamically accessible 
perturbations. We prove that equilibria without signature changes are struc- 
turally stable and those with changes are structurally unstable. In Sec. [6] we 
define critical states of the linearized Vlasov equation that are structurally 
unstable under perturbations that are further restricted. We prove that a 
mode with the opposite signature of the continuum is structurally unstable 
and that the opposite combination cannot exist unless the system is already 
unstable. Finally, in Sec. [7], we conclude. 



2 Noncanonical Hamiltonian form 



The class of equations of interest have a single dependent variable ((x,v,t), 
such that for each time t, £ : T> — > R, where the particle phase space T> is some 
two-dimensional domain with coordinates (x,v). The dynamics is assumed 
to be Hamiltonian in terms of a noncanonical Poisson bracket of the form 



{F,G} 



dxdv ( 



v 



5F_ SG 



(1) 



where [/, g] := f x g v — f v g x is the usual Poisson bracket, where the subscripts 
denote partial differentiation, and 5F/5( denotes the functional derivative of 
a functional F[Q. The equation of motion is generated from a Hamiltonian 
functional H[C] as follows: 



<Zt = {C,n} = -[(,e\, 



(2) 



where £ := 5T-L/5C,. The Poisson bracket of ([!]) is noncanonical: it uses only 
a single noncanonical variable (, instead of the usual canonically conjugate 
pair; it possesses degeneracy reflected in the existence of Casimir invariants, 
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C = f v dxdvC((), that satisfy {F, C} = for all functional F; but it does 
satisfy the Lie-algebraic properties of usual Poisson brackets. For further 
details see £21 123 HD1 EU ■ 

For the Vlasov-Poisson equation we assume T> — X x R, where X C R 
or X = S, the circle, the distinction will not be important. The dependent 
variable is the particle phase space density f(x,v,t) and the Hamiltonian is 
given by 

nn = \ [ dx[dw 2 f + l [ dxi^w (3) 

z J x Jr 1 J x 

where is a shorthand for the functional dependence on / obtained through 
solution of Poisson's equation, <f) xx = 1 — Lfdv, for a positive charge species 
with a neutralizing background. Using 8%/S f = £ = v 2 /2 + 0, we obtain 

ft = {f, H} = - [f, 8} = -vf x + <p x f v , (4) 

where, as usual, the plasma frequency and Debye length have been used to 
nondimensionalize all variables. 

This Hamiltonian form for the Vlasov-Poisson equation was first pub- 
lished in |16j . For a discussion of a general class of systems with this Hamil- 
tonian form, to which the ideas of the present analysis can be applied, see 
[21]. In a sequence of papers [2U [lj [TH [251 121 [26] various ramifications of the 
Hamiltonian form have been explored - notably, canonization and diagonal- 
ization of the linear dynamics to which we now turn. 

Because of the noncanonical form, linearization requires expansion of the 
Poisson bracket as well as the Hamiltonian. Equilibria, (o, are obtained by 
extremization of a free energy functional, F = H + C, as was first done for 
Vlasov-like equilibria in [27]. Writing C = Co + Ci an d expanding gives the 
Hamiltonian form for the linear dynamics 

( lt = {( u H L } L (5) 

where the linear Hamiltoinian, Hl = \ f v dxdv CiO £i, is the second variation 
of F, a quadratic from in d defined by the symmetric operator O, and 
{F, G}l = f v dxdv ( [F 1 ,Gi\ with F\ := 5F/5C,\. Thus the linear dynamics 
is governed by the time evolution operator T ■:=—{■ , H L } L = [£ , O ■ ]. 

Linearizing the Vlasov-Poisson equation about an homogeneous equilib- 
rium, fo(v), gives rise to the system, 

fit = -vfi x + hxfo ( 6 ) 

(t>ixx = ~ dvfi, (7) 
Jr 
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for the unknown f\(x, v, t). Here f := dfo/dv. This is an infinite-dimensional 
linear Hamiltonian system generated by the Hamiltonian functional: 

H L [h} = -\ [ dx [ dv ^ \j\\ 2 + I [ dx\<p lx \ 2 . (8) 
A J x Jr Jo a J x 

We concentrate on systems where x is an ignorable coordinate, and either 
Fourier expand or transform. For Vlasov-Poisson this gives the system 

i f r 

f kt = -ikvf k + — ^ dv fk{v, t) =: -T k f k , (9) 
ft Jr 

where fk(v,t) is the Fourier dual to fi(x,v,t). Perturbation of the spec- 
turm of the operator defined by Eq. ^ is the primary subject of this paper. 
The operator is a Hamiltonian operator generated by the Hamiltonian 
functional 

H L [f k , U] = \ (- f dv j \fk\ 2 + \<P k \ 2 ) , (10) 

u \ JR. J / 



k 

with the Poisson bracket 



Observe from (11) that k E N and thus fk and f-k are independent variables 
that are almost canonically conjugate. Thus the complete system is 

fkt = -T k fk and f- kt = -T-uf-k , (12) 
from which we conclude the spectrum is Hamiltonian. 

Lemma 2.1 If X is an eigenvalue of the Vlasov equation linearized about 
the equilibrium fv{v), then so are —A and A (complex conjugate). Thus if 
A = 7 + iu, then eigenvalues occur in the pairs, ±7 and ±iu, for purely 
real and imaginary cases, respectively, or quartets, A = ±7±iw, for complex 
eigenvalues. 



Proof That —A is an eigenvalue follows immediately from the symmetry 
T_fc = — Tk, and that A is an eigenvalue follows from Tkfk = —(Tkfk)- D 
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In [U [TSJ [21 ESI ^ was shown how to scale /& and to make them 
canonically conjugate variables. In order to do this requires the following 
definition of dynamically accessibility, a terminology introduced in [28, 29J. 

Definition A particle phase space function k is dynamically accessible from 
a particle phase space function h, if g is an area-preserving rearrangement 
of h; i.e., in coordinates k(x,v) = h(X(x,v),V(x,v)), where [X, V) = 1. A 
peturbation 5h is linearly dynamically accessible from h if 5h = [G, h], where 
G is the infinitesimal generator of the canonical transformation (x, v) <H- 
(X,V). 

Dynamically accessible perturbations come about by perturbing the par- 
ticle orbits under the action of some Hamiltonian. Since electrostatic charged 
particle dynamics is Hamiltonian, one can make the case that these are the 
only perturbations allowable within the confines of Vlasov-Poisson theory. 

Given an equilibrium state /o, linear dynamically accessible perturbations 
away from this equilibrium state satisfy Sfo = [G, fo] = G x /q. Therefore 
assuming the initial condition for the linear dynamics is linearly dynamically 
accessible, we can define 

qk(v,t) = f k and Pk{v,t) = -if- k /(kf ) (13) 

without worrying about a singularity at the zeros of f' Q and k = 0. With the 



definitions of (13), the Poisson bracket of (11) achieves canonical from 



{F,G} L = }^ / dv — — — . (14) 

f-^ Jr V oQk 0p k opt dq k J 

The full system has the new Hamiltonian % = % + UV in frame mov- 
ing with speed U, where V = j v dxdvvf . Linearizing in this frame yields 
the linear Hamiltonian = Hi + P^, from which we identify the linear 
momentum 

i r t. 

P L [f k J- k ] = - 9 Yl dv-\f k \\ (15) 
z k=1 Jr Jo 

which must be conserved by the linear dynamics. It is easy to show directly 
that this is the case. 



Lemma 2.2 The momentum Pl defined by (15) is a constant of motion, 
i.e., {P L ,H L } = 0. 
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Proof This follows immediately from (12): j R dv {fkT-k + f-kTk) = 0. □ 



Observe, that like the Hamiltonian, H^, the momentum Pi is conserved 
for each k, which in all respects appears only as a parameter in our system. 
Assuming the system size to be L yields k = 2irn/L with n e N, and, 
thus, this parameter can be taken to be in M + /{0}. Alternatively, we could 
suppose X — R, Fourier transform, and split the Fourier integral to obtain an 



expression similar to (11) with the sum replaced by an integral over positive 
values of k. For the present analysis we will not be concerned with issues 
of convergence for reconstructing the spatial variation of fi(x, v, t), but only 
consider k e IR + /{0} to be a parameter in our operator. We will see in 
Sec.[3]that the operator possesses a continuous component to its spectrum. 
But, we emphasize that this continuous spectrum of interest arises form the 
multiplicative nature of the velocity operator, i.e. the term vfk of T^, not from 
having an infinite spatial domain, as is the case for free particle or scattering 
states in quantum mechanics. In the remainder of the paper, / will refer 
to either fi or which will be clear from context, and the dependence on 
k will be suppressed, e.g. in T^, unless k dependence is being specifically 
addressed. 



3 Spectral stability 

Now we consider properties of the evolution operator T defined by rt9j). We 
define spectral stability in general terms, record some properties of T, and 
describe the tools necessary to characterize the spectrum of T. We suppose 
fk varies as exp(— iut), where u is the frequency and ico is the eigenvalue. 
For convenience we also use u := u/k, where recall k e M + . The system 
is spectrally stable if the spectrum of T is less than or equal to zero or 
the frequency is always in the closed lower half plane. Since the system is 
Hamiltonian, the question of stability reduces to deciding if the spectrum is 
confined to the imaginary axis. 

Definition The linearized dynamics of a Hamiltonian system around some 
equilibrium solution, with the phase space of solutions in some Banach space 
£>, is spectrally stable if the spectrum a(T) of the time evolution operator T 
is purely imaginary. 

Spectral stability does not guarantee that the system is stable, or that the 
equilibrium f is linearly stable. (See e.g. [ID] for general discussion). The 
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solutions of a spectrally stable system are guaranteed to grow at most sub- 
exponentially and one can construct a spectrally stable system with polyno- 
mial temporal growth for certain initial conditions. (See e.g. [30J for analysis 
of the Vlasov system. ) 

Spectral stability relies on functional analysis for its definition, since the 
spectrum of the operator T may depend on the choice of function space B. 
The time evolution operators arising from the types of noncanonical Hamil- 
tonian systems that are of interest here generally contain a continuous spec- 
trum [21] and the effects of perturbations that we study can be categorized 
by properties of the continuous spectrum of these operators. In general for 
the operators of [2T], the operator T is the sum of a multiplication operator 
and an integral operator. In the Vlasov case, the multiplicative operator 
is if - and the integral operator is f' Q Jdv-. As we will see, the multiplica- 
tion operator causes the continuous spectrum to be composed of the entire 
imaginary axis except possibly for some discrete points. 

Instability comes from the point spectrum. In particular, the linearized 
Vlasov Poisson equation is not spectrally stable when the time evolution 
operator has a spectrum that includes a point away from the imaginary axis, 
with the necessary counterparts implied by Lemma 2|2 For the operator 



T this will always be a discrete mode; i.e. an eigenmode associated with an 
eigenvalue in the point spectrum. 

Theorem 3.1 The one-dimensional linearized Vlasov-Poisson system with 
homogeneous equilibrium fo is spectrally unstable if for some k e M + and u 
in the upper half plane, the plasma dispersion relation 

e(k, u) := 1 - AT 2 t dv^- = . 

Jr v ~ u 

Otherwise it is spectrally stable. 

Proof The details of this proof are given in plasma textbooks. It follows 
directly from ^ and ([7]), and the assumption fi ~ exp(ikx — iut). □ 

Using the Nyquist method that relies on the argument principle of com- 
plex analysis, Penrose [31] was able to relate the vanishing of e(k, u) to the 
winding number of the closed curve determined by the real and imaginary 
parts of s as u runs along the real axis. Such closed curves are called Penrose 
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plots. The crucial quantity is the integral part of e as u approaches the real 
axis from above: 



lim — 

u^0+ IT 




H[f ]{u)-if Q {u), 



v — u 



where H[f'^\ denotes the Hilbert transform, H[f'^\ = -J-dv fb/(v — u), where 
j := PV L indicates the Cauchy principle value. (See [32] for an in depth 
treatment of Hilbert transforms.) The graph of the real line under this map- 
ping is the essence of the Penrose plot, and so we will refer to these closed 
curves as Penrose plots as well. When necessary to avoid ambiguity we will 
refer to the former as e-plots. 

For example, Fig. [T] shows the derivative of the distribution function, 
/o, for the case of a Maxwellian distribution and Fig. [2] shows the contour 
H\fo\ — ifo(u) that emerges from the origin in the complex plane at u = — oo, 
descends, and then wraps around to return to the origin at u = oo. From 
this figure it is evident that the winding number of the e(k, w)-plot is zero 
for any fixed fcel, and as a result there are no unstable modes. 

Making use of the argument principle as described above, Penrose ob- 
tained the following criterion: 

Theorem 3.2 The linearized Vlasov-Poisson system with homogeneous equi- 
librium /o is spectrally unstable if there exists a point u such that 



with /q traversing zero at u. Otherwise it is spectrally stable. 

Penrose plots can be used to visually determine spectral stability. As 

2 

described above, the Maxwellian distribution /o = e~ v is stable, as the 
resulting e-plot does not encircle the origin. However, it is not difficult 
to construct unstable distribution functions. The superposition of two dis- 
placed Maxwellian distributions, f = e~^ +c ^ + e~^~ c ^ , is such a case. As 
c increases the distribution goes from stable to unstable. Figures [3] and [4] 
demonstrate how the transition from stability to instability is manifested in 
a Penrose plot. The two examples are c = 3/4 and c = 1. (Note, the nor- 
malization of /o only affects the overall scale of the Penrose plots and so is 
ignored for convenience.) It is evident from Fig. [4] that for some fcsK the 
£-plot (which is a displacement of the curve shown by multiplying by —k~ 2 




and 




u 



>0, 
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Figure 1: f for a Maxwellian distribution. 




Figure 2: Stable Penrose plot for a Maxwellian distribution. 
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-1.5 -1 -0.5 0.5 1 1.5 2 2.5 3 



mm 

Figure 3: Penrose plot for a stable superposition of Maxwellian distributions 




-1.5 -1 -0.5 0.5 1 1.5 2 2.5 3 



H[fS\ 

Figure 4: The unstable Penrose plot corresponding to two separated Maxwell 
distributions. 
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and adding unity) will encircle the origin, and thus will be unstable for such 
k- values. 

We now are positioned to completely determine the spectrum. For con- 
venience we set k — 1 when it does not affect the essence of our arguments, 
and consider the operator T: f \-> ivf — if^J f in the space W 1 ' 1 ^), but 
we also discuss the space L X (]R). The space H^ 1 ' 1 (IR) is the Sobolev space 
containing the closure of functions under the norm ||/||i,i = ||/||i + H/'Hi- 
Thus it contains all functions that are in L 1 (IR) whose weak derivatives are 
also in L 1 (M). First we establish the expected facts that T is densely defined 
and closed. 

In W 1,1 the operator T is the sum of the multiplication operator and a 
bounded operator - that it is densely defined and closed follows from the 
fact that the multiplication operator is densely defined and closed in these 
spaces, where 

D 1 (T) := {f\vf G W 1 ' 1 (R)}. 

Theorem 3.3 The operator T : W 1 ' 1 (R) ->■ W^QR) with domain £>i( T ) is 
both (i) densely defined and (ii) closable. 

Proof (i) The set of all smooth functions with compact support, C£°(R) is 
a subset of D\. This set is dense in W /1 ' 1 (M) so D\ is dense and T is densely 
defined, (ii) The operator T is closable if the operator v is closable because T 
and v differ by a bounded operator. The multiplication operator v is closed 
if for each sequence f n C W 1,:L (R) that converges to either vf n converges 
to or vf n does not converge. Suppose vf n converges. At each point f n 
converges to 0. Therefore vf n converges to at each point, so vf n converges 
to if it converges. □ 

Therefore there exists some domain D such that the graph (D,TD) is 
closed. 

In determining the spectrum of the operator T, denoted cr(T), we split 
the spectrum into point, residual, and continuous components as follows: 

Definition For A G (t(T) the resolvent of T is R(T, A) = (T — A/) -1 , where 
/ is the identity operator. We say A is (i) in the point spectrum, cr p (T), if 
T — XI fails to be injective, (ii) in the residual spectrum, a r (T), if R(T, A) is 
not densely defined, and (iii) in the continuous spectrum, <r c (T), if R(T,\) 
is densely defined but unbounded. 
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Using this definition we characterize the spectrum of the operator T. 



Theorem 3.4 The component cr p (T) consists of all points A = iu G C 
where 1 — k~ 2 J R dv f^/iy — u) = 0, cr c (T) consists of all A = iu with u G 
icr p (T) flK), and cr r (T) contains all points A = iu in the complement of 
c p (T) U cr c (T) that satisfy /o( n ) = 0- 

Proof By the Penrose criterion we can identify all the points in the point 
spectrum. If 1 — k~ 2 L dvf Q /{v — u) = then iu = A G c p (T). Because 
the system is Hamiltonian these modes will occur for the linearized Vlasov- 
Poisson system in quartets (two for and two for TL&), as follows from 
Lemma 2][2 It is possible for there to be discrete modes with real frequencies 



and these will occur in pairs. If for real u the map u \- > e passes through the 
origin then there will be such an embedded mode. 

For convenience we drop the wavenumber subscript k on /& and add 
the subscript n to identify /„ as an element of a sequence of functions that 
converges to zero with, for each n, support contained in an interval of length 
2e(n) surrounding the point u and zero average value. Let u G M and choose 
the sequence {/„} so that e(n) — > 0. Then for each n 



\\R(T,iu)\\ > 



n 1,1 



> 



(? - U) /n|| 1,1 

II /"II 1,1 



\v — ttll^i.l^-e^+ejll/nlll.l 
1 



\\V — U\\w^(u-e,u+e) 

In the above expression, W l,1 {u — e,u + e) refers to the integral of |/| + |/'| 
over the interval (u — e,u + e). Therefore the resolvent is an unbounded 
operator and iu — A is in the spectrum. If the frequency u has an imaginary 
component ry then iu)\\ < I/7 so unless iu — A is part of the point 

spectrum it is part of the resolvent set. 

The residual spectrum of T is contained in the point spectrum of T*. 
The dual of W 1 ' 1 is the space W~ l)1 defined by pairs (g, h) G W~ 1,x with 
HO, ||_i,i < 00 (cf. p3j). The operator T*(g,h) = i(vg - h + J\gfo - 
hfH )dv , — vh) is the adjoint of T. If we search for a member iu = A of the 
point spectrum we get two equations, one of which is (v — u)h = 0. This 
forces h = because h cannot be a 5-function in W~ 1,1 . The other equation 
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is then (v — u)g + J gf' Q dv = which can only be true if the integral is zero 
or if (v — u)g is a constant. For this g = and the resulting equation for 
u is the same equation as that for the frequency of the point modes of T . If 
the integral is zero then g = 5(v — u) is a solution when /q(u) = 0. Therefore 
the residual spectrum contains the points A = in satisfying f^iu) — 0. □ 

This characterization of the spectrum fails in Banach spaces with less 
regularity than W" 1 ' 1 , such as L p spaces, because the Dirac 5 is not con- 
tained in the dual space. In this case the residual spectrum vanishes because 
cr p (T*) = cr p (T). This calculation is nearly identical to that of Degond [30J, 
who characterizes the residual spectrum slightly differently than we do. In 
any event, the result is that the Penrose criterion determines whether T is 
spectrally stable or not. If the winding number of the e-plot is positive, then 
there is spectral instability and if it is zero there is spectral stability. 

4 Structural stability 

Spectral stability characterizes the linear dynamics of a nonlinear Hamilto- 
nian system in a neighborhood of an equilibrium. The main question now is 
to determine when a spectrally stable system can be made spectrally unstable 
with a small perturbation. When this is impossible for our choice of allowed 
perturbations, we say the equilibrium is structurally stable, and when there 
is an infinitesimal perturbation that makes the system spectrally unstable 
we say that the equilibrium is spectrally unstable. We can make this more 
precise by stating it in terms of operators on a Banach space. 

Definition Consider an equilibrium solution of a Hamiltonian system and 
the corresponding time evolution operator T for the linearized dynamics, 
with a phase space some Banach space B. Suppose that T is spectrally 
stable. Consider perturbations 5T of T and define a norm on the space of 
such perturbations. Then we say that the equilibrium is structurally stable 
under this norm if there is some 5 > such that for every \\8T\\ < 5 the 
operator T + ST is spectrally stable. Otherwise the system is structurally 
unstable. 

Because we are dealing with physical systems it makes sense to have 
physical motivation for the choice of norm on the space of perturbations. In 
this paper we are interested in perturbations of the Vlasov equation through 
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changes in the equilibrium. This choice is motivated by the Hamiltonian 
structure of the equations and Krein's theorem for finite-dimensional systems. 
In general the space of possible perturbations is quite large, but perturbations 
of equilibria give rise to operators in certain Banach spaces and motivate the 
definition of norm. Even in the case of unbounded perturbations there may 
exist such a norm (see Kato [B], for instance). 

Consider a stable equilibrium function f . We will consider perturbations 
of the equilibrium function and the resulting perturbation of the time evo- 
lution operator. Suppose that the time evolution operator of the perturbed 
system is T + ST. In the function space that we will consider these pertur- 
bations are bounded operators and their size can be measured by the norm 
||5T||. This norm will be proportional to the norm of ||5/o||, where Sfo is the 
perturbation of the equilibrium. 

Definition Consider the formulation of the linearized Vlasov-Poisson equa- 
tion in the Banach space W /1 ' 1 (IR) with a spectrally stable homogeneous equi- 
librium function f . Let Tf 0+ sf be the time evolution operator corresponding 
to the linearized dynamics around the distribution function fo + Sfo. If there 
exists some S depending only on f such that Tf 0+ sf is spectrally stable when- 
ever ||5Ta/ || = \\Tf — Tf 0+ gf \\ < S, then the equilibrium f is structurally 
stable under perturbations of f . 

The aim of this work is to characterize the structural stability of the 
linearized Vlasov-Poisson equation. We will prove that if the perturbation 
function is some homogeneous Sfo and the norm is W 1 ' 1 (and L 1 as a con- 
sequence) every equilibrium distribution function is structurally unstable to 
an infinitesimal perturbation in this space. This fact will force us to consider 
more restricted sets of perturbations. 

4.1 Winding number 

We need to compute the winding number of Penrose plots and the change in 
winding number under a perturbation, both in this section and the rest of the 
paper. We use the fact that one way to compute the winding number is to 
draw a ray from the origin to infinity and to count the number of intersections 
with the contour accounting for orientation. 

Lemma 4.1 Consider an equilibrium distribution function f . The winding 
number of the Penrose e-plot around the origin is equal to s S n {f'o ( u )) f or 
all u G M~ , satisfying f (u) = 0. 
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To calculate the winding number of the Penrose e-plot using this lemma 
one counts the number of zeros of /q on the negative real line and adds them 
with a positive sign if /q is positive, a Penrose crossing from the upper half 
plane to the lower half plane, a negative sign if /q is negative, a crossing from 
the lower half plane to the upper half plane, and zero if u is not a crossing 
of the x-axis, a tangency. This lemma comes from the following equivalent 
characterization of the winding number from differential topology pi]. 

Definition If X is a compact, oriented, /-dimensional manifold and / : X — > 
M. l+1 is a smooth map, the winding number of / around any point z G — 
f(X) is the degree of the direction map u: X — > S l given by u(x) = jyjfpfj- 

In our case the compact manifold is the real line plus the point at oo 
and 1 = 1. The degree of u is the intersection number of u with any point 
on the circle taken with a plus sign if the differential preserves orientation 
and a minus sign if it reverses it. The lemma is just a specialization of 
this definition to the negative x-direction on the circle. If more than one 
derivative of ft vanishes at a zero of ft there is a standard procedure for 
calculating the winding number by determining if there is a sign change in 
ft at the zero. 



4.2 Structural instability of general /q 

In a large class of function spaces it is possible to create infinitesimal pertur- 
bations that make any equilibrium distribution function unstable. This can 
happen in any space where the Hilbert transform is an unbounded operator. 
In these spaces there will be an infinitesimal 8ft such that H[8ft] is order 
one at a zero of f' Q . Such a perturbation can turn any point where ft = 
into a point where H[ft + 8 ft] > as well. Because 8ft is small and the 
region where H[Sft] is not small is also small, the only effect on the Penrose 
plot will be to move the location of the zero. Thus, such a perturbation will 
increase the winding number and cause instability. 

We will explicitly demonstrate this for the Banach space H /1 ' 1 (IR) and, by 
extension, the Banach space L 1 n Cq. This will imply that any distribution 
function is infinitesimally close to instability when the problem is set in 
one of these spaces, implying the structural instability of every distribution 
function. 

Suppose we perturb ft by a function 6 ft. The resulting perturbation to 
the operator T is the operator mapping / to 8ft J dvf. In the space W 1,1 this 
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is a bounded operator and thus we take the norm of the perturbing operator 
to be || 5/q || 1,1. Now we introduce a class of perturbations that can be made 
infinitesimal, but have Hilbert transform of order unity. 
Consider the function x( v , h, d, e) defined by 



X 



hv/e 
h sgn(-u) 

h + d/2 + e/2-v/2 
-h - d/2 - e/2-v/2 




\v\ < e 

e < \v\ < d + e 
2h + d + e>v>d + e 
2h + d + e>-v>d + e 
\v\ > 2h + d + e 



Figures [5] and [6] show the graph of x an d its Hilbert transform, respec- 
tively. 
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Figure 5: The perturbation x for e = e 10 , h = d = .1. 



Lemma 4.2 If we choose d = h and e = e^^^ h \ then for any 5, 7 > 
we can choose an h such that ||xl|i,i < ^ an d fdvx/ v > 1 — O(h), and 
\j-dv x/{u — v)\ < \-f/u\ for \u\ > \2h + d+ e\. 

Proof In the space W 1,1 the function x nas norm 2h 2 + 2hd + he + 4h, 
which is less than any S for small enough h. We can compute the value of 
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Figure 6: The Hilbert transform of X- 



the Hilbert transform of this function at a given point u by calculating the 
principal values: 



fdv 



X 



v — u 



hu f\u — e\ 
e \ \u + e 



+ h\og 



\d + 


e — u\ 


\u + d + e 




e + u\ 


|e — it 





+ -(d + e + 2/i- w)lo£ 



\d + e + - u\ 

\d + e — u\ 
\d + e + 2h + u\ 



+ l{d + e + 2h + u)log 

2 \ \d + e + u 



(16) 



We analyze the asymptotics of this function as h, d, and e go to zero, with the 
desiderata that i) the norm of x goes to zero, ii) the maximum of the Hilbert 
transform of x is 0(1), and iii) there is a band of vanishing width around the 
origin outside of which the Hilbert transform can be made arbitrarily close 
to zero. 



Note that (16) can be written as a linear combination of translates of the 
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function x log x: 



dv 



X 



((« - e) logflw - e|) - (u + e) log(|w + e|)) 



h 

e 

-(d + u + e) logfld + u + e|) - -(d - u + e) logfld -u + e\) 

1 



+ -(d + u + e + 2h)\og(\d + u + e + 2h\) 
+ \(d - u + t + 2h)\og(\d - u + e + 2h\) . 



(17) 



The function xlogx has a local minimum for positive x at x = 1/e. This 
is the point at which the function is most negative. It has zeros at x = 
and sc = 1. For values of it, d, e, /i close to zero all of the arguments of the 
log functions are less than 1/e. Therefore, for \u\ < d + e + 2h the xlogx 
terms are all monotonically decreasing functions of the argument x. Of the 
terms of (17), -((« — e) log(|w — e|) — (u + e) log(|w + e|)) has by far the 
largest coefficient as long as e is much smaller than h. We choose h = d 
and e = 0(e _1 / /l ). Then the terms that do not involve e are all smaller than 
(6h + e) log(6/t + e). With these choices x satisfies 

x (0) = 2-(h + e- 1 / h )log(\h + e- 1/h \) + (3h + e- 1 / h )log(\3h + e- 1/h \) 
= 2 + 0(h\ogh). 

Consider the pair of functions — (u + c) log(|w + c|) + (u — c) log(|w — c|). The 
derivative with respect to u is — log(|u + c|) + log(|u — c|). This is zero for 
w = and for u > it is always negative and the pair is always decreasing, 
and for small values of h the pair is guaranteed to be positive. Suppose that 
u > e. Then we can bound the term with the h/e coefficient: 



h 



\{u-e) log(|u - e|) -(u + e) logflu + e| 



h 



log 



\u 



h 



< 



[u - e) log 



u + e 
1 - 



1 + - 

u 

-|(w-e) log(e- e/ " 
e 

h(u — e) 



— 2elog(|ti 

— 2elog(|-u - 



+ 2|/ilog(|u + e|)| 



u 



+ 2 \h\og(\u 
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For u » e, for example if u = 0(h 2 ), this term is O(hlogh). Therefore, for 
\u\ > h 2 we have x = O(hlogh) which can be made arbitrarily small. When 
\u\ > 3h + e the function \ decreases at least as fast as 0(l/u). With these 
choices of h, d, and e, the norm of \ is 0(h), which proves the Lemma. □ 

Now we state the theorem that any equilibrium is strucutrally unstable 
in both the spaces W 1,1 and L 1 fl C . 

Theorem 4.3 A stable equilibrium distribution fo G C 2 is structurally un- 
stable under perturbations of the equilibrium in the Banach spaces W 1 ' and 

Proof If f Q is stable then the Penrose e-plot of /q has a winding number 
of zero. Because the point at oo corresponds to a crossing where f' Q goes 
from negative to positive there exists a point Uq with Jq(uq) — that is an 
isolated zero, iff/oK^o) < 0, and fo(u ) < 0. Let F = sup |/o'l- Choose h 
to always be smaller than the distance from uq to the nearest of f' Q . Then 
if e = 0(e~ 1 / /l ) and d = h the support of x( u ~~ u o) will contain only one 
zero of /q. For h small enough the slope of x a ^ u o w ih be greater than F 
so that the function f' Q + \ wm be positive for u in the set (uq, u + ) for some 
u + in the support of x- Similarly J'q + x wiU be negative for u in the set 
(u~,uq) for some u~ in the support of x- Because x has compact support 
the function f' Q + x is positive in a neighborhood outside of the support of x 
so that the intermediate value theorem guarantees one additional zero of the 
function f' Q + x for u > u and also for u < u . Choose x so that this Hilbert 
transform of f' Q + x is positive at the point uq and h small enough that it is 
negative before the next zero of /q + x on either side of uq. Then the winding 
number of f + x is positive because an additional positive crossing has been 
added on the negative real line. 

Because the norm of x is 0(h) in both W 1 ' 1 and L l the distribution f 
is unstable to an arbitrarily small perturbation and is therefore structurally 
unstable. □ 

Thus we emphasize that we can always construct a perturbation that 
makes our linearized Vlasov-Poisson system unstable. For the special case 
of the Maxwellian distribution, Fig. [7] shows the perturbed derivative of the 
distribution function and Fig. [8] shows the Penrose plot of the unstable per- 
turbed system. Observe the two crossings created by the perturbation on the 
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positive axis as well as the negative crossing arising from the unboundedness 
of the perturbation. 



In some sense Theorem |4.3| represents a failure of our class of perturba- 
tions to produce any interesting structure for the Vlasov equation. Indeed 
signature appears to play no role in delineating bifurcation to instability. In 
order to derive a nontrivial result we develop a new theory analogous to the 
finite-dimensional Hamiltonian perturbation theory developed by Krein and 
Moser. This new theory involves a restriction to dynamically accessible per- 
turbations of the equilibrium state. This is natural since the noncanonical 
Hamiltonian structure can be viewed as the union of canonical Hamilto- 
nian motions (on symplectic leaves) labelled by the equilibrium state - to 
compare with traditional finite-dimensional theory requires restriction to the 
given canonical Hamiltonian motion under consideration. 
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Figure 7: f Q + x f° r a Maxwellian distribution. 



5 The Krein-Moser theorem 

For linear finite-dimensional Hamiltonian systems, Hamilton's equations are 
a set of first order linear ODEs. If the Hamiltonian is time-independent, then 
the behavior of solutions is characterized by the eigenfrequencies. If all of the 
eigenfrequencies are on the real axis and nondegenerate, then the system will 
be stable. If there are degenerate eigenvalues the system will be stable as long 
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Figure 8: Penrose plot for perturbed Maxwellian. 

as the time evolution operator does not have any nontrivial Jordan blocks, 
but there will be secular growth if it does. Any complex eigenfrequencies will 
lead to instability The Hamiltonian of a linear finite-dimensional Hamilto- 
nian system is a quadratic form in the canonical variables. If we consider 
perturbations of the coefficients of the quadratic form it is trivial to define a 
notion of small perturbations, as the resulting perturbation of the Hamilto- 
nian will be a bounded operator. Krein and Moser independently proved a 
theorem characterizing the structural stability of these systems in terms of 
a signature, a quantity that amounts to the sign of the energy evaluated on 
the eigenvector of a mode. We consider a special case of their theorem. 

Theorem 5.1 (Krein-Moser) Let H define a stable linear finite- dimensional 
Hamiltonian system. Then H is structurally stable if all the eigenfrequencies 
are nondegenerate. If there are any degeneracies, H is structurally stable if 
the assosciated eigenmodes have energy of the same sign. Otherwise H is 
structurally unstable. 

This Krein-Moser theorem gives a clear picture of the behavior of these 
systems under small perturbations. The eigenfrequencies move around, but 
remain confined to the real line unless there is a collision between a positive 
energy and negative energy mode, in which case they may leave the axis. This 
theorem was first proved by Krein in the early 1950s and later rediscovered 
by Moser in the late 1950s. Our goal is to place the perturbation theory 
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of infinite-dimensional Hamiltonian systems in the language of the finite- 
dimensional theory. 

The appropriate definition of signature for the continuous spectrum of 
the Vlasov-Poisson equation was introduced in [TJ [2] (see also [26J), where 
an integral transform was also introduced for constructing a canonical trans- 
formation to action-angle variables for the infinite-dimensional system. The 
transformation is a generalization of the Hilbert transform and it can be 
used to show that the linearized Vlasov-Poisson equation is equivalent to the 
system with the following Hamiltonian functional: 

oo „ 

Hl = ^2 dua k (u)u k (u)J k (u,t) , (18) 
k=i 

where Uk(u) = \ku\ and cr fe (u) = — sgn(fcu/g(u)) is the analog of the Krein 
signature corresponding to the mode labeled by u G K. (Note, the transfor- 
mation can always be carried out in a frame where /q(0) = 0. Because the 
Hamiltonian does not transform as a scalar for frame shifts, which are time 
dependent transformations, signature is frame dependent. The Hamiltonian 
in a shifted frame is obtained by adding a constant times the momentum 



of (15) to Hi. Later we will see that Hamiltonians that can be made sign 



definite is some frame are structurally stable in a sense to be defined.) 

Definition 1 Suppose /q(0) = 0. Then the signature of the point u e R is 
-sgn(u/o(u)). 

Below is a graph that illustrates the signature for a bi-Maxwellian distri- 
bution function. 



5.1 Dynamical accessibility and structural stability 

Now we discuss the effect of restricting to dynamically accessible pertur- 
bations on the structural stability of /q. In this work we only study per- 
turbations of fo that preserve homogeneity. Because dynamically accessible 
perturbations are area preserving rearrangements of fo, it is impossible to 
construct a dynamically accessible perturbation for the Vlasov equation in a 
finite spatial domain that preserves homogeneity. 

To see this we write a rearrangement as (x,v) -H- (X,V), where V is a 
function of v alone. Because [X, V] = 1 and V(v) is not a function of x, we 
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Figure 9: Signature for a bi-Maxwellian distribution function. 

have V'dX/dx = 1, or X = x/V. If the spatial domain is finite, this map 
is not an diffeomorphism unless V — 1. In the infinite spatial domain case, 
this is not a problem and these rearrangments exist. First we note that a 
rearrangement cannot change the critical points of /q. 

Lemma 5.2 Let (X, V) be an area preserving diffeomorphism, and let V be 
homogeneous. Then the critical points of f (V) are the points l /_1 (w c ) ; where 
v c is a critical point of fo(v). 

Proof By the chain rule df (V(v))/dv = V{v)' f^(V(v)). The function V ^ 
because (X, V) must be a diffeomorphism. Therefore the critical points 
occur when /q(V) = or at points v = ^ _1 (t> c ). □ 

Consider the perturbation \ that was constructed earlier. If v c is a nonde- 
generate critical point of /o such that fo(v c ) < 0, then we want to prove that 
there is a rearrangement V such that fo(V) = fo(v) + f_ x( v ' ~ v c )dv' or 
that dfo(V)/dv = fb(v) + x(v — v c ). Such a rearrangement can be constructed 
as long as the parameters defining x, the numbers h, d, e, are chosen such that 
fo(v) +x( v ~ v c) has the same critical points as fo(v). Using Morse theory it 
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is possible to construct a V so that fo(V) = fo(v) + / x + 0((v — v c ) 3 ), where 
0((v— f c ) 3 ) has compact support and is smaller than fo(v) — fd (v c )(v— v c ) 2 /2. 

Theorem 5.3 Let v c be a nondegenerate critical point of fo with fb(v c ) < 0. 
Then there exists a rearrangement V such that fo(V) = fo(v) + f x + 0((v — 
v c ) 3 , where O is defined as above. 

We omit the proof but it is a simple application of the Morse lemma. In 
order to apply the Morse lemma fo must be C 2 . This is not restrictive for 
practical applications where typically fo is smooth. The rearrangement of fo 
can also be made to be smooth if desired. 

Using this result we prove a Krein-like theorem for dynamically accessible 
perturbations in the W 1,1 norm. 

Theorem 5.4 Let fo be a stable equilibrium distribution function for the 
Vlasov equation on an infinite spatial domain. Then f is structurally stable 
under dynamically accessible perturbations in W 1 ' 1 , if there is only one solu- 
tion of fb(v) = 0. If there are multiple solutions, f is structurally unstable 
and the unstable modes come from the zeros of f' Q that satisfy fd(v) < 0. 

Proof Suppose that f has only one zero on the real line. Because f is 
an equilibrium this zero will have /q > 0. Because a dynamically accessi- 
ble perturbation can never increase the number of critical points, it will be 
impossible to change the winding number of the Penrose plot to a positive 
number. Therefore fo is structurally stable. 

Suppose that fb = has more than one solution on the real axis. Using 
the previous theorem perturb f' Q by x( v ~ v c) m a neighborhood of a critical 
point v c with fo(v c ) < 0. This will increase the winding number to 1 since it 
will add a positively oriented crossing on the negative real axis for the correct 
choices of h, d, and e in the definition of x- The norm of x can be made as 
small as necessary and therefore f is structurally unstable. Since no other 
critical points with f ' < can be created the only critical points that lead 
to instabilities are the ones that already exist having f ' < 0. □ 

The implication of this result is that in a Banach space where the Hilbert 
transform is an unbounded operator the dynamical accessibility condition 
makes it so that a change in the Krein signature of the continuous spectrum 
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is a necessary and sufficient condition for structural instability. The bifurca- 
tions do not occur at all points where the signature changes, however. Only 
those that represent valleys of the distribution can give birth to unstable 
modes. 



6 Krein bifurcations in the Vlasov equation 

We identify two critical states for the Penrose plots that correspond to the 
transition to instability. In these states the system may be structurally un- 
stable under infinitesimal perturbations of f' Q in the C n norm for all n. The 
first critical state corresponds to the existence of an embedded mode in the 
continuous spectrum. If the equilibrium is stable, then such an embedded 
mode corresponds to a tangency of the Penrose plot to the real axis at the 
origin. If the system is perturbed so that the tangency becomes a pair of 
transverse intersections, then the winding number of the Penrose plot would 
jump to 1 and the system would be unstable. Considering a parametrized 
small perturbation, we see that the value of k for the unstable mode, will 
correspond to some value of k ^ for which the embedded mode exists. 



Figures 10 and 11 illustrate a critical Penrose plot for a bifurcation at k ^ 0. 



We explore this bifurcation in Sec. 6.1 



Another critical state occurs when H[f'^\ = at a point where f trans- 
versely intersects the real axis. If the Hilbert transform of f is perturbed, 
there will be a crossing with a negative H[Jq], and the winding number will 
be positive for some k. This mode enters through k = because the smaller 
the perturbation of H[f'^\ the smaller k must be for to be unstable. Figure 



12 is a critical Penrose plot corresponding to the bi-Maxwellian distribution 
with the maximum stable separation. We explore this kind of bifurcation in 
Sec. 16.31 



6.1 Bifurcation at k ^ 

The linearized Vlasov equation can support neutral plasma modes embedded 
within the continuous spectrum. The condition for existence of a point mode 
is the vanishing of the plasma dispersion relation on the real axis, 

e(u) = l-H[ti}+i& = 0. (19) 
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Figure 12: Critical Penrose plot for a bi-Maxwellian distribution function. 

If the spatial domain is unbounded the point modes will be analogues of the 
momentum eigenstate solutions of the Schrodinger equation and have infinite 
energy. Any violation of the Penrose criterion will guarantee the existence of 
zeros of the plasma dispersion function on the real axis because k can take 
any value in this case. 

If the plasma dispersion vanishes at some u and f[l(u) = 0, there is an 
embedded mode in the continuous spectrum. The signature of the continuous 
spectrum will not change signs at the frequency of the mode and we will 
extend the definition of signature to the point u even though f^iu) = 0. The 
signature of an embedded mode is given by sgn(udeR/du) (see [H[25]). The 
signature of the continuous spectrum is — sgn(tt/o). These signatures are the 
same if the value of f' Q in a neighborhood of its zero is the same sign as H\f^\. 

We will prove that if /o is stable and mildly regular, it is impossible 
for there to be a discrete mode embedded in the continuous spectrum with 
signature that is the same as the signature of the continuous spectrum sur- 
rounding it. The proof has a simple conceptual outline. Suppose that there 
exists a discrete mode with the same signature as the continuum. Then there 
exists some point u satisfying f' {u) = 0, — sgn(/g) = sgn(<9e#/<9w) in a neigh- 
borhood of u, and H[f'^\{u) = 0. Perturbations of f centered around this 
point will give the Penrose plot a negative winding number, contradicting 
the analyticity of the plasma dispersion function in the upper half plane. We 
need f to be Holder continuous so that the Penrose plot is continuous and 



28 



for the plasma dispersion function to converge uniformly to its values on the 
real line. 

Lemma 6.1 Let g be a function defined on the real line such that g is Holder 
and let h = H[f]. Then the functions g z , h z that are the solutions of the 
Laplace equation in the upper half plane satisfying f z = f and g z = g on the 
real line converge uniformly to f and g . 

Proof Because g can be defined as a bounded and continuous function on 
R U {00} and the g z are analytic, the g z must converge uniformly to g. The 
same properties hold for h and h z must converge to h — H[g\. □ 

Lemma 6.2 Let f be the derivative of an equilibrium distribution function 
and let f' Q be sufficiently regular such that the assumptions of the previous 
lemma are true. Then the Penrose plot that is associated with f Q cannot 
have a negative winding number. 

Proof The Penrose plot associated with f Q is the image of the real line under 
the map e(u) = 1 — H[fo\ + ifb- This is naturally defined as an analytic 
function if u is in the upper half plane. By the argument principle the image 
of R + it under this map has a non-negative winding number. Both the real 
and imaginary parts of this map converge uniformly to their values on the 
real line. Therefore the Penrose plot is a homotopy of these contours, making 
it possible to parametrize the contours by some t such that the distance from 
the Penrose plot to the contour produced by the image of R + i5 is always less 
than some rj(t) that goes to 0. If the winding number of the Penrose plot were 
negative, there would be some t for which the winding number was negative 
because the winding is a stable property under homotopy, contradicting the 
analyticity of the map. □ 

Theorem 6.3 Let f' Q and /q be Holder continuous. If f is stable there are 
no discrete modes with signature the same as the signature of the continuum. 

Proof Because f is stable the winding number of the Penrose plot is equal 
to 0. Assume that there is a discrete mode with the same signature as 
the contiuum surrounding it. Then there exists a point u with /q(m) = 0, 
fo(u) = 0, and sgn(fl ) (u+5)dH[f'^\(u)/du) = 1. Then we search for a function 
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g such that the Penrose plot of /q + g has a negative winding number. If 



such a function exists it will contradict Lemma 6.1 Because /q is Holder 
den/du is bounded away from zero in a neighborhood of the point. Suppose 
that in this neighborhood there is only one zero of f' Q . Then define g such 
that g has one sign, is smooth and has compact support, and such that 
the \dH[g]/du\ < |/q| in this neighborhood. Then for small enough g the 
function f^ + g will have two zeros in a neighborhood of the point. Then 
both of the crossings will correspond to crossings of negative orientation and 
the resulting winding number will be —1, a contradiction. □ 

Corollary 6.4 If f is stable it is impossible for there to be a point where 
f' Q = 0, f» < 0, and H[Q > 0. 

If fo is unstable the winding number is positive. In this case it may be 
possible for modes with the same signature as the continuum to exist. It is 
possible for a positive energy mode to be embedded in a section of negative 
signature and a negative energy mode to be embedded in a section of positive 
signature. This situation is structurally unstable under perturbations that 
are bounded by the C n norm and remains so even when a linear dynamical 
accessibility constraint is enforced. 

Theorem 6.5 Let f Q be the derivative of an equilibrium distribution function 
with a discrete mode embedded in the continuous spectrum. Then there exists 
an infinitesimal function with compact support in the C n norm for each n 
such that /q + 5f is unstable. 

Proof Suppose that H[fH] is non-zero in a neighborhood of the embedded 
mode. Define a dynamically accessible perturbation 5f = hf . Then assume 
that /q" 7^ at the mode. If we define h such that it does not vanish at the 
mode we find that Sf" = h" f' Q + h'f({ + hf{j' and therefore we can choose h 
such that the discrete mode becomes a crossing. This can be done with h 
infinitesimal and smooth. The resulting perturbation will have an infinites- 
imal effect on f' Q . The new crossings will cause a violation of the Penrose 
criterion, and therefore the system with the embedded mode is structurally 
unstable. □ 

This is an analog of Krein's theorem for the Vlasov equation for the case 
where there is a discrete mode. As a result of this we see that all discrete 
modes are either unstable or structurally unstable. 
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6.2 Little-big man theorem 



Consider a linearized equilibrium that supports three discrete modes. The 
signature of each mode depends on the reference frame. There is a result 
that applies to a number of Hamiltonian systems, the three-wave problem in 
particular [USE]), that gives a condition on the signature of the modes and 
their frequency in some reference frame such that no frame shift can cause 
all the modes to have the same signature. In a shifted frame the Hamiltonian 
changes by virtue of the frequencies in the action-angle from being doppler 
shifted. Sometimes such shifts can render the Hamiltonian sign definite. The 
little-big man theorem for finite systems indicates when this cannot happen. 
Such a result is also possible for the point spectrum of the Vlasov equation. 

Theorem 6.6 Let f' Q be the derivative of an equilibrium distribution func- 
tion that has three discrete modes (elements of the point spectrum) with real 
frequencies. Consider a reference frame where all of the modes have posi- 
tive frequency. Then represent the energies of the three modes as a triplet 
(± ± ±) where the plus and minus signs correspond to the signature of each 
mode, with the first mode being the one with the lowest frequency and the last 

the one with the highest frequency. Then if the triplet is of the form (H h) 

or ( — I — ) there is no reference frame in which all the modes have the same 
signature. If the triplet has any other form, then there is a reference frame 
in which all the modes have the same signature. 

The formula for the energy of an embedded mode is sgn(ude R /du). If 
there are three embedded modes in a frame where the frequencies are all 
positive the triplet is 
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Proof If this is (H h) then as we shift frames the possible triplets are 

(0-+), ( h), (-0+), (- + +), (- + 0), (- + -). All of these are indefinite. 

The other possibile initially indefinite triplet is ( h). However if we shift 

the two — modes to negative frequency the triplet becomes (+ + +). All 
other examples are either definite or reduce to one of these two. □ 
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6.3 Bifurcation at k = 



Assume that there are no embedded modes and that fo is stable, but that 
there is a point that has fo = and H[f \ = 0. This is the critical state for a 
bifurcation at k — 0. This can be destabilized in the same way as the critical 
state for k ^ 0. There will be a perturbation that makes H[f' Q ] < without 
changing /q at that point. Therefore the Penrose plot becomes unstable and 
the equilibrium is structurally unstable. 

Theorem 6.7 Suppose that f' Q is a stable equilibrium distribution function 
that has a zero at u of both f and H[f'^\. Then f is structurally unstable 
under perturbations bounded by the C n norm for all n. 

Proof Let 5h be symmetric about the point u, be smooth with compact 
support and have its first n derivatives less than some e. Then let 5f' Q = 
— H[5h]. The resulting perturbation to H[f' Q ] is h. If h is positive at u, then 
by the symmetry of h f' Q + 5f' Q has a zero at u and H[f ] + h is positive 
there. Thus the Penrose plot has a positive winding number and is unstable. 
Therefore f is structurally unstable. □ 

The previous two sections demonstrated that when the Penrose plot is 
critical, no amount of regularity is sufficient to prevent fo from being struc- 
turally unstable. However, when the Penrose plot is not critical all that is 
required is that a small perturbation only change the Penrose plot by a small 
amount in addition to a condition to prevent perturbations near v — oo. 
Suppose we arbitrarily restrict the support of the perturbations so that 
\v\ < v max . Then if we increase the required regularity such that sup(//"[5/o]) 
is bounded there will be some 5 such that for all Sf^ with ||5/o|| < 5 the 
distribution f + Sfo is structurally stable. This restriction can be moti- 
vated physically by restricting the particles in the distribution function to be 
travelling slower than the speed of light. 

7 Conclusion 

We have considered perturbations of the linearized Vlasov-Poisson equation 
through changes in the equilibrium function. The effect of these perturba- 
tions on the spectral stability of the equations is determined by the class of 
allowable perturbations and the signature of the contiuous and point spectra. 
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Every equilibrium can be made unstable by adding an arbitrarily small func- 
tion from the space W 1,1 . If we rearrange fo then only when the signature of 
fo changes sign can an arbitrarily small perturbation destabilize it. When fo 
is stable discrete modes always have the opposite signature of the spectrum 
surrounding them. The equilibria are structurally unstable under C n small 
perturbations for all n. The signature of the spectrum and the signature of 
the discrete mode can never be the same. 

This generalization of Krein's theorem is more complicated than the 
finite-dimensional original. However the basic ideas of Krein's theorem are 
still important in the infinite-dimensional case. When the perturbations are 
more restricted than just belonging to W 1 ' 1 the structural stability is deter- 
mined by the signature of the spectrum. Just as in Krein's theorem there 
must be a positive signature interacting with a negative signature to produce 
structural instability. 

This paper was devoted primarily to the Vlasov equation, but other non- 
canonical Hamiltonian systems admit to a similar treatment, e.g. the 2D 
Euler equation with shear flow equilibria, and we hope to chronicle such 
cases in future publications. 
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